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The interplay of electron-phonon (el-ph) and electron-electron (el-el) interactions in epitaxial 
graphene is studied by directly probing its electronic structure. We found a strong coupling of 
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electrons to the soft part of the Ai„ phonon evident by a kink at 150±15 meV, while the coupling 
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of electrons to another expected phonon E2 9 at 195 meV can only be barely detected. The possible 
role of the el-el interaction to account for the enhanced coupling of electrons to the Ai g phonon, 
and the contribution of el-ph interaction to the linear imaginary part of the self energy at high 



binding energy are also discussed. Our results reveal the dominant role of the A± g phonon in the 
el-ph interaction in graphene, and highlight the important interplay of el-el and el-ph interactions 
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in the self energy of graphene. 
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FIG. 1: (Color online) (a) ARPES data taken along KM direction (solid line in the inset), (b) 
Dispersion (black curve) extracted by fitting the raw data. The dashed line is the fit using two 
straight lines with different slopes. Within 20 meV below ~Ep, the dispersion is affected by the 
resolution, and therefore we fit the dispersion only in the range between -250 meV and -20 meV. The 
gray dotted line is a guide for the deviation of the low energy dispersion from the extrapolation of 
the high energy dispersion, (c) Extracted velocity as a function of energy from dispersions plotted 
in Panel b. The symbols are the raw data and the black solid line is a guide for the eye. (d) 
Extracted real part of the self energy as a function of energy Re£(E)=E(k)-v&kp where E(k) is the 
measured dispersion and kp is the Fermi wave vector. The inset shows the MDC width. 

Electron-phonon (el-ph) coupling is among the most important interactions, since it is at 
the origin of a variety of interesting phenomena, such as the hopping-like charge transport in 
organic semiconductors [TJ [2] , charge density wave formation [3] , metal-insulator transition, 
superconductivity [I] and ballistic transport [SJ E|. The el-ph interaction is particularly 
intriguing in graphitic materials, where the special electronic properties of Dirac fermions 
and the interplay of el-ph and electron-electron (el-el) interaction result in a wide range of 
novel physics [7j. Because of its peculiar point-like Fermi surface, which can be connected by 
the wave vectors of the phonons at T and K [8] , electron screening of the lattice vibrations 
decreases dramatically around these two points, causing two Kohn anomalies [5]. Moreover 
in the case of single layer graphene, the peculiar band structure also results in the breakdown 
of the Born Oppenheimer approximation [TUl ITT] ; a shift of the E2 9 phonon frequency as a 
function of carrier concentration [TOl [TT] IT2| IT3] and sample thickness [Hj; and a predicted 
anomalous phonon-induced self energy [151 EH H3 that deviates from that of conventional 
metals [18J. 
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Despite the intense research effort, two key components in understanding el-ph interac- 
tion - the phonon modes involved and the coupling strength of the interaction, have not 
been resolved. Angle-resolved photoemission spectroscopy (ARPES) is an ideal tool in this 
respect, as it directly measures the renormalized electronic band structure of a material and 
therefore provides direct insights about many-body interactions. In recent years ARPES has 
been successfully used to detect the signature of the el-ph interaction in the electronic spec- 
tra in the form of a kink, in both graphite [191 E3 IZE] and graphene [22~1 [23] . However, not 
only is there a discrepancy in the value of the observed coupling strength [TjJl 12DI [2D [221 [23] 
with respect to the theoretical predictions (HE EE EZ], but also consensus on which and 
how many phonon modes are involved has been missing so far. Theoretically it was pro- 
posed that due to the Kohn anomalies at T and K, both the F,2 g (195 meV) and Ai s (165 
meV) phonons contribute to the el-ph interaction [15] . Experimentally although a kink has 
been reported in the electronic dispersion [221 123] , the large uncertainty in the kink energy 
makes it difficult to distinguish which, if not both, phonons are involved. Therefore a more 
detailed study with improved data quality is needed to complete our understanding of the 
el-ph coupling in graphene and provide key insights for the el-ph coupling in other graphitic 
materials. 

In this letter we present a high resolution ARPES study of the el-ph interaction and its 
contribution to the electron self energy in epitaxial graphene. The greatly improved data 
quality with reduced noise level has enabled us for the first time to nail down the kink energy 
in the electronic dispersion to 150±15 meV and to reveal additional fine structures in the 
electron self energy at 60±15 meV and 200±15 meV. More importantly, the direct compar- 
ison between the electronic dispersion measured here and the reported phonon dispersion 
relation [8] has allowed us to identify the soft part of the Ai g phonon (Kohn anomaly) as 
the main scattering channel responsible for the ARPES kink, and the fine structure at 200 
meV in the self energy with the E 2g mode. The enhanced coupling to the A ls mode with 
respect to the E2 9 mode together with the much larger experimental el-ph coupling strength 
A~0.14 as compared to theoretical one is discussed in terms of Coulomb interactions. In 
addition, we report the linear imaginary part of the self energy at high binding energy with 
similar magnitude along various directions, which reflects the contribution from both el-ph 
and el-el interactions. Our results point to the dominant role of A l9 phonon in the el-ph 
interaction in graphene, and highlight the important interplay of el-ph and el-el interaction 
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FIG. 2: (color online) (a) Intensity maps for the neighboring K and K 7 points as a function of 
energy and momentum. The two arrows show the el-ph interaction with the phonon and K\ g 
phonon respectively, (b) Phonon dispersions for the Ai s near K and E2 9 near r [U L2| . 

in the intriguing physics of Dirac fermions in graphene. 

High resolution ARPES data were taken on single layer epitaxial graphene at Beamline 
12.0.1 (Figs. 1-3) and Beamline 7.0.1 (Fig. 4) of the Advanced Light Source (ALS) of the 
Lawrence Berkeley National Laboratory with a total energy resolution of 25 meV and 35 
meV respectively. Samples were grown on a n-type SiC wafers as detailed elsewhere [2U 122] • 
The samples were measured with 50 eV photon energy at a temperature of 25K, with vacuum 
better than 3.0 xl0~ n Torr. 

Figure la shows an ARPES intensity map taken through the Dirac point (K point). One 
can easily identify a characteristic energy (pointed to by a horizontal black arrow), where the 
slope of the dispersion (i.e. velocity) changes and the intensity suddenly decreases due to the 
disappearance of coherent peaks in the energy distribution curves (EDCs) at high binding 
energy. These are typical signatures of electron-boson coupling, where the broadening of 
the spectra beyond the kink energy is due to the onset of the bosonic mode self energy 
[T9| |20~1 |2T| [22| 123] . To identify the exact bosic modes involved in the coupling and the 
strength of the coupling, we extract the dispersion relation from the peak positions and the 
ImS from the peak width by fitting the momentum distribution curves (MDCs). The high 
statistics of the data in panel b allows us to nail down the kink position to 150±15 meV [26J. 
The extracted kink energy is also consistent with a sudden drop of the MDC width (inset 
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FIG. 3: (Color online) (a) Dispersions along various directions (labeled as a, (3, 7 and 8) shown 
in the Brillouin zone in the inset. The dashed lines are fits of the dispersion using two straight 
lines with different slopes. The dotted line is a guide to the eye for the deviation at low binding 
energy. The horizonal gray shadow highlights the kink energy, (b) Corresponding ReS extracted 
by subtracting the bare band dispersion from the measured dispersion. The bare band dispersion 
is taken as a straight line connecting the dispersions at and -250 meV. 

of panel d) and a change of the electron velocity (panel c), both occurring at -150 meV. 
From the renormalization of the electron velocity we extract the el-ph coupling constant A 
= Vft/vjr-l, where v&=1.0x 10 6 m/s is the bare band velocity, and vp the renormalized Fermi 
velocity vp=0.87x 10 6 m/s. This gives an experimental A = 0.14 which is almost an order of 
magnitude larger than the predicted value of 0.02 for the A\ g phonon [T3j. The identification 
of the kink at 150 meV with strength of 0.14 is also supported by recent scanning tunneling 
microscope measurements |27j . 

To check whether other phonon modes contribute to this large value of A, we show in 
panel d the real part of the electron self energy, ReS. In addition to the main peak at -150 
meV that dominates ReS, two additional fine structures at pa -60 meV and ~ -200 meV can 
also be resolved. The existence of these fine structures indicates the involvement of other 
collective modes in the coupling. Since the area underneath the ReS is an indication of 
the coupling strength, clearly these additional modes contribute only a small fraction to the 
total coupling constant, and hence cannot be responsible for the large discrepancy between 
the experimental and theoretical A. 

To single out the allowed scattering processes, in figure 2 we compare the ARPES constant 
energy map at Ep and at the kink energy (panel a) with the predicted phonon dispersion 
(panel b) [8] . Clearly the soft part of the K\ g phonon near the zone corner K point is the only 
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FIG. 4: (Color online) (a-c) Dispersions along various directions. The inset shows the constant 
energy map at -3 eV, where the trigonal distortion can be clearly observed. The three lines in the 
inset label the three cuts shown in panels a-c. (d-f) Extracted MDC width as a function of energy 
for data shown in panels a-c. The gray line is a guide for the transition from a more linear behavior 
to a more quadratic behavior, (g-i) Imaginary part of the self energy Im£=^v-Ak, where v is the 
velocity at each energy and Ak is the MDC width subtracted by the width at E^? (dotted lines in 
panels d-f) to take care of the impurity scattering and finite experimental resolution. 



mode with the right energy and momentum q^|ri^| to scatter states separated in energy 
by 150 meV (kink energy) from near the K point to the K' point (inter- valley scattering), 
therefore being likely the dominant source for the kink in the dispersion and the maximum 
peak in the ReS. Similarly, the E2 5 phonon near the T point has the right energy (195 meV) 
and momentum q^0 to connect states between 200 meV and the Fermi energy within the 
same K point (intra-valley scattering) and is responsible for the fine structure in ReS at 
200 meV. Since the sample is slightly electron doped, the q vector for intra- and inter-valley 
scattering is ~4%|rif| larger than q=|r/f | and q=0, but still in the proximity of the two 
Kohn anomalies [9j. Finally, the fine structure in ReS at 60 meV is likely due to coupling 
with an out-of-plane phonon as reported by STM studies [28] . 

Figure 3 compares the el-ph interaction along different directions. A similar kink is 
present in the dispersion along all the directions at the same energy of 150 ± 15 meV (see 
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gray region in panel a). Although similar additional fine structures involving the two other 
phonon modes are also observed in the self energy in panel b, the most important finding is 
that the A ls phonon is still the dominant one. 

Although on a qualitative level the data presented here are in good agreement with theo- 
retical calculations, on a quantitative level there are two key differences. 1) Experimentally 
we found that the intervalley scattering with the A\ g phonon is by far the dominant scatter- 
ing source. This is in contrast to theoretical prediction where both intervalley (A lg phonon) 
and intravalley (E 2g phonon) scatterings are treated in an almost equal footing, although the 
latter is decreased by half [PJ. 2) The experimental el-ph coupling strength of 0.14, mostly 
accounted for by the Ai s phonon (as discussed in figure 1), is much larger than the predicted 
value of 0.02 [15]. This holds even if finite experimental resolution, which makes A twice 
as big, is taken into account [15]. Therefore, additional mechanism needs to be included 
to explain the observed enhancement, by approximately a factor of 3, of the el-ph coupling 
strength. One likely candidate is through the interplay with el-el interaction, as pointed 
out theoretically [T51 UH] - More specifically, it has been argued that this interplay gives rise, 
in the presence of a linear dispersion, to a linear imaginary part of the self energy ImE in 
agreement with experimental reports [22j [29] . and that the el-ph coupling contributes to 
1/3 of its total magnitude [16]. By extending this study to the entire momentum region 
(Figure 4), we shown that this linearity in ImE survives with a similar magnitude through- 
out the entire Dirac cone, even when the dispersion is not linear because of the trigonal 
distortions (Fig.4(b,c)), and is hence a general property of Dirac fermions. Panels (a-c) 
show the ARPES intensity maps along three different directions. From TK (panel a) to 
MK direction (panel c), both the extracted dispersion (see solid black line in panels a-c) 
and MDC width between -1 and -3 eV change from a linear to a quadratic behavior, due 
to the trigonal distortion (see also inset of panel a) [30]. The corresponding ImE are shown 
in panels (g-i). Clearly, even when the dispersion is not linear (panel c), ImE still shows a 
linear dependence with similar magnitude (panel i), suggesting that the overall contribution 
of the el-ph interaction to the self energy of graphene is comparable along all directions, in 
line with the isotropic el-ph coupling reported in Fig. 3 and theoretical prediction [31]. 

These data clearly establish the importance of the interpaly between el-ph and el-el 
interaction suggesting that the latter might be responsible for the observed enhancement of 
the coupling strength. Indeed it has been argued that in the completely unscreened case, the 
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el-el interaction can enhnace the coupling to the K\ g phonon near K by up to a factor of 3, 
leaving the coupling to the E 2g phonon near T almost unaffected [32] ■ This picture reconciles 
the disagreement between the theoretical and experimental coupling strength measured by 
ARPES and can also account for the large intensity ratio between the 2D and 2D' peaks 
reported by Raman [33], which is likely due to the enhanced renormalization of the el-ph 
coupling of peak D. Finally it is interesting to note that a similar enhanced coupling to 
phonons with non-zero wave vectors through el-el interaction also occurs in the case of 
transition metal dichalcogenides where the quasiparticles are Dirac fermions, resulting as 
well in a linear ImS [31]. These similarities suggest that the physics here discussed is not 
only a property of graphene but a more general property of Dirac materials. 

In conclusion, we have reported the strong interplay of el-ph and el-el interactions in 
graphene. We identified the dominant role of the Ai g phonons at 150 meV in the el-ph 
interaction along the various directions near the K point. Although the fine structures due 
to coupling with other phonon modes are observed, we show that the enhancement of the 
coupling to the soft part of the A lg mode is likely induced by the interplay between el-el 
and el-ph interactions. This study demonstrates the important role of this interplay in a 
Dirac fermion system, and highlights the importance of including both interactions in the 
self energy of graphene. 
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